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We extend a previous analysis of spatial correlation functions for classical electromagnetic vector 
fields near a perfectly conducting boundary [PRE, 73, 036604 (2006)] to the case of an isotropic 
semi-infinite medium with planar interface and characterized by a first-order impedance boundary 
condition. The analytical results are illustrated with calculations for the case of point separations in 
the direction perpendicular to the interface. For the incident plus reflected fleld, the dependence of 
the complex-valued and inhomogeneous spatial correlation function on the permittivity, permeability 
and conductivity of the medium is determined. For the refracted field, the spatial correlation is again 
complex-valued but homogeneous and highly sensitive to the value of the refractive index. Based 
on the derived dependencies, novel nonlocal measurement methods for precision characterization of 
electromagnetic material properties are suggested. The influence of the directionality of incidence 
for electromagnetic beams is investigated. Narrowing the beam width results in a slower decrease 
of the amplitude of the correlation function as a function of point separation. Previously obtained 
asymptotic results for statistically homogeneous random free fields are retrieved as special cases. 
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I. INTRODUCTION 



By extending a series of earlier studies for homoge- 
neous free fields we recently derived and analyzed 
spatial correlation functions of inhomogeneous random 
classical electromagnetic (EM) fields in . In that anal- 
ysis, the configuration consisted of a perfect electrically 
conducting (PEC) infinite planar boundary, resulting in 
a superposition of incident and reflected hemispherical 
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statistically isotropic random fields in front of the inter- 
face. Since a PEC boundary exhibits constant, i.e., angle- 
and polarization-independent reflection coefficients for 
both perpendicular (transverse electric, TE) and parallel 
(transverse magnetic, TM) wave polarizations, the reflec- 
tion of a statistically isotropic incident fleld exhibits an 
isotropic angular spectrum as well. On the other hand, 
the spatial correlation of the incident plus reflected flelds 
was found to be inhomogeneous (i.e., dependent on the 
absolute distance of either one of the two point locations 
with respect to the interface), as a consequence of the 
statistical field anisotropy imposed by the EM boundary 
condition at the interface. First-order statistics (prob- 
ability distributions) of the energy density for this con- 
figuration were derived in which were also found to 
exhibit inhomogeneity through action at-a-distance. 

In this paper, we extend this previous study by con- 
sidering spatial correlation functions for EM fields in the 
presence of a semi-infinite isotropic medium, as a sec- 
ond canonical configuration of fundamental interest. The 
impedance boundary condition causes the reflection and 
transmission coefficients to depend on both the polariza- 
tion state (TE, TM, or hybrid) and the angle of incidence. 
As a result, the angular spectra of both the reflected 
and refracted flelds are now no longer hemispherically 
isotropic, but are nonuniformly weighted across the solid 
angle of incidence. Secondly, we investigate the influence 
of directional incidence (sectorial solid angle of incidence 
centered around a central direction), including narrow 
EM beams as a limiting case. 

The present analysis and results are relevant to sev- 
eral practical problems of interest, e.g., coherence prop- 
erties of stellar light transmitted through an atmosphere 
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or inside optical instruments, radio waves reflected by the 
Earth's soil or ionosphere, multipath scattering by man- 
made objects or precipitation, multi-mode cavities, etc. 
The field coherency - which is the basic EM quantity in 
such scenarios - is expressed in terms of reflection and 
transmission coefficients for plane waves impinging onto 
a single planar interface. By extension, results more gen- 
eral multi-layer configurations are obtained without diffi- 
culty, by simply substituting the Fresnel coefficients for a 
semi-infinite medium with corresponding expressions for 
stratified media. 



II. THEORY 



A. Reflected plus incident fields 




1^, s zi-ir j 



Consider a semi-infinite isotropic medium with scalar 
permittivity e — er^o, permeability fi — /irAio, conduc- 
tivity a and first-order surface impedance rj = -^////e = 

770 \/ [i^rl^r where |arg(7/)| < 7r/4. This medium occu- 
pies the half-space z < (Fig. [T|). We assume a time- 
harmonic random incident field (E*, H*) which can be ex- 
panded as an isotropic angular spectrum of plane waves, 
each specified by a triplet (f*,H*,k*) and propagating 
toward the interface, i.e., k* • 1^ < 0, where ||k*|| = 
fco — LUy/JIoeo is the (constant) free-space wavenumber of 
each incident plane wave. A harmonic time dependence 
exp{jLut) is assumed and suppressed. The overall inci- 
dent electric field at r = r* can then be represented 
as 



E* (r^) = — 



£\n) exp (-jk' • r*) dfi. 



(1) 



The integral ^ is valid for inhomogeneous random fields 
and hence applicable to the present configuration, un- 
like homogeneous random fields which, strictly, require 
a Fourier-Stieltjes representation incorporating general- 
ized fields (distributions). Incidence and refraction of 
the plane waves in the upper and lower hemispheres Qq 
and f2 is governed by angles respectively, for 

their propagation direction relative to the surface nor- 
mal. The wavenumber within the refracting medium is 
I Ik* 1 1 = k = uj^ITlI ~ kn^/Urer. The refracted electric 



field E(r) is expanded in a similar way as ([T]), mutatis 
mutandis. 

For a general stratified multi-layered medium, in- 
cluding the particular case of a semi-infinite isotropic 
medium, TE and TM waves constitute a set of uncoupled 
eigenmodes. Hence, their contributions to the resultant 
field can be evaluated separately and then superimposed, 
at any location. We refer to Sections II and III of ^] for 
notations and detailed calculations of the TE/TM de- 
composition of a random field with respect to the surface 
normal. For an isotropic semi-infinite medium and TE 
polarization, the resultant (i.e., incident plus refiected) 
field at ro(a;o, j/Oi -Zq) within the local plane of incidence 



FIG. 1: Coordinate system and local plane of incidence (<^o ~ 
0, l</,o = ly) for single TE wave component reflected and 
refracted by a semi-inflnite isotropic medium. 



= 00 is 

£y cxp(-jk • ro) = £ya cxp (-jfcofi'o sin 6*0) 

X {[1 + r_L(6'o)]cos(fcozocos6'o) 
-r^(6lo)]sin(fco2;o cos^o)} 

(2) 

£ 

'Hxexp{-jk- ro) = — cose^o exp (-j/coPo sin6'o) 

X { [1 - F^ (6*0)] cos (fco^^o cos 6*0) 
+j [1 + r^iOo)] sin (fco^o cos 6*0)} 

(3) 

£ 

Hz exp(-jk • ro) = — sin 6*0 exp (-jfeoPo sin 6*0) 

X {[1 + Fj^(6'o)]cos(fco2:o cos 6*0) 
-r^(6'o)]sin(fcoZo cos 6*0)} 

(4) 

where go = a;ocos0o + 2/0 sin 0o, £yo = ^^^cos^o + 
fffg cos^o sin^o, = 'S^oCOsV'o, = -^osin^o, in 
which the random polarization angle ■00 is uniformly dis- 
tributed within the local transverse plane spanned by 
l0o and Igg. In ©-Q, k = k* when k • 1^ < and 
k = k*^ = (I - 21zlz) ■ k' for k • 1^ > 0. The Fresnel TE 
reflection coeflacient is 



r.(^o) = ^fe"-^oVf;-^^ + ^^^, (5) 



rjku + rjQ ^fc2 — /cq + k^u'^ 

in which u = cos^Oj where = denotes a deflnition. Com- 
paring ©-gl) above with (4)-(5) in 0, it follows that the 
z-dependence of the resultant fleld is no longer spatially 
harmonic when 7^ ^±{do) 7^ ±1, unlike in the case of 
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a PEC surface. Similarly, for the wave components that 
are TM with respect to the plane of incidence, we have 



where 



exp(-jk • ro) = 



Sxex.p^-jk- ro) 



£:^exp(-jk- ro) = 



-Hyo exp 

x{[i-r 
+j [1 + r 

rioHy„ cos 

x{[i + r 
+j [1 - r 

rjoHyg sin ( 

x{[i-r 
+.1 [1 + r 



i-jkogo sin 6*0) 
(6*0)] cos(fcozocos6'o) 
(6*0)] sin(fco2;o cos 6*0)} 
(6) 

?oexp (-j/cofiio sin 6*0) 
(6*0)] cos(fcozocos6'o) 
(6*0)] sin(fcoZo cos 6*0)} 
(7) 

'0 exp(-jfco£'o sin 6*0) 
(6*0)] cos(fco2ocos6'o) 

(0o)] sin(fcoZo cos 6*0)} 
(8) 



lyl 



ly2 



lyS - 



lyi - 



2C f\l-T^{u)\' 
Jo 

X sin (fco^iw) sin (fcoZ2w) (H) 

2C f li + r^Hi' 

Jo 

X cos{kQZiu) cos{kQZ2u) du (12) 

j2C f\i-rAu)] [i + riH] 

Jo 

X sm.{koZiu) COS {koZ2u) du (13) 

-j2C f\l + T^{u)] [i-rKu)] 
Jo 

X COS (fcoziw) sin (fco2;2w) du (14) 



with Tiy 



Ti-d 



cos 6*0 sin 00, Ti. 



00 



TYosinV'o, Weo = HqCostAo, 
TM reflection coefficient 



Sq/tiq, and Fresnel 



%) = 



Tjy'k^ 



fcg 



(9) 



Whilst we shall limit the further analysis to point separa- 
tions in normal direction (Arl^ = Aziz), the results are 
easily extended to arbitrary directions using the method- 
ology outhned in Section IV of 

For the TE waves, substitution of ^ into ([T]), eval- 
uated at two locations ri,2 = 21,212 for ro, enables 
the calculation of Ey {ri) ■ E*{r2) via double integration 
with respect to corresponding ranges f2i and f22 [bI, eq. 
(16)], where the asterisk denotes complex conjugation. 
This is followed by ensemble averaging of this product, 
assuming delta-correlated random field components @, 
eq. (17)], i.e., {£i{ni) ■ £*{n2)) = {Siei^i) ■ £*2g{^2)) + 
{£i^{ni) ■ £^^2)) - 2CS[{ni U ^2) \ {fli n ^2)], where 

C = (|£oP)/4- If, in addition, each complex Carte- 
sian component exhibits a zero mean value, then these 
impositions on the first- and second-order moments de- 
fine, unambiguously, a 3-D complex (6-D real) multi- 
variate Gauss normal distribution with independent and 
identically distributed components The results also 
apply to more general distributions for £1^2, provided 
that their first- and second-order moments satisfy the 
stated expressions. It follows that for a general isotropic 
impedance boundary condition, the TE field coherency 
{Ey{zi)E*{z2)) can be written as a sum of four terms. 



Note that /ys 7^ I*^ and 9(/j/3 + lyA) 7^ unless zi =22, 
so that the spatial coherencies are in general complex- 
valued. The integrals pT|) - P^ evaluate to closed-form 
but cumbersome expressions. Expressions for the TM co- 
herencies follow in an analogous manner by substituting 
and dH]) into ([I]), yielding 



{E^{zi)E*{z2)) = /xl + 1x2 + 1x3 + 1x4 



where 



Ixi = 2C / |l-r||(M)| 
Jo 

xu^ sin (koziu) sin {koZ2u) du 
pi 

1x2 = 2C / \l + Ti\{uf 
Jo 

xu^ COS {koziu) cos {koZ2u) du 

1x3 = j2C f [i-r||(w)] [i-f r|(u) 

Jo '- 
xu^ sin (koZiu) cos {koZ2u) du 

1x4 = -32c [ [i + r||(u)] h-rUu 

Jo 

xu COS (fco^iw) sin (fco2;2M) du 



(15) 



(16) 



(17) 



(18) 



(19) 



and 



{Ey{zi)E*{z2)) = lyl + Iy2 + Iy3 + Iy4 



(10) 



{Ezizi)E;{z2)) = hi + 1,2 + Iz3 + h4 



(20) 
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where 



hi = 2C / i + r||(u) 



X (l — u^) sin (fco^iu) sin (fcoZ2M) du (21) 
pi 

Iz2 = 2C7 / |l-r||(u)|' 
Jo 

X (l — u^) cos (fcoziu) cos {koZ2u) du (22) 

j2C /'[i + r||(7.)] [i-q(«)" 

X (l — u ) sin (fco^iw) cos (fcoZ2u) dw (23) 

-j2C / [i-r||(u)] [i + r|(u) 

X (l — u^) cos (/cqZiu) sin {koZ2u) du (24) 

For the normal field = E^lz, the tangential field Ej = 
Exlx + Eyly, and the total, i.e., vector field E = E^lx + 
Eyly + Ezlz, we have 



Iz4 = 



{Ez{zi)El{z2)) = 
{Et{z{)E;(z2)) 
{E{zi)E*{z2)) 



a=x,y e=l 

E E^- 

a— a:,y,z 



(25) 
(26) 
(27) 



B. Refracted fields 

For the field transmitted (refracted) across the inter- 
face, we obtain in an analogous manner, with the aid of 
the field transmission coefficients Tl || (6*0) and Snell's law 
ko sin 6q ~ k sin 6, 




{Ey{zi)E;{z2))=2C / \T^{u)\^ 
Jo 



X exp 



{Ex{z{)e:{z2))^2c \mu)Y 

Jo 



du (28) 




X exp ]kn[zi - Z2)\l I — ) - 1 + du (29) 

{EM)E:iz2)) = 2c£ miu^ [jj (1 - u^) 
jkoizi ~ Z2)\ ( ^] -l + u2 du (30) 



X exp 



where the Fresnel TE and TM transmission coefficients 
are 



T^{u) = 1 + T^{u) 



2rjku 



rjku + 770 a/A:^ — /cq + k^u- 



= , (31) 



cos (^0 
COS 9 



:i + r||(u)] 

2riku 



klu'^ 



(32) 



respectively. The coherencies now exhibit complex- 
harmonic dependencies on the separation distance only, 
as in the case of free random fields in an infinite homo- 
geneous medium. Thus, unlike for the incident plus re- 
ficcted field, the spatial correlation of the refracted field 
is homogeneous ^ i.e., dependent on fcAz — k\z\ — Z2\ only. 
Physically, this is a consequence of the fact that no in- 
terference exists beyond the interface. Nevertheless, be- 
cause of the ^o-dependence of Tiji, the coherency of the 
refracted field is different from that of the incident field. 



III. NUMERICAL RESULTS 

A. Good conductor 

For good but imperfect nonmagnetic conductors (cr ^ 
tjeg, e — eoi A* = Mo), we can approximate 77/770 ~ 
^ uje^l {2a)(\ ± j) and ~ 0. The refiection coefficients 
then become 



V^(u) 



yu - 770 

r]u + r]o' 



rii(^) 



V - Vou 

77 + 77o7i 



(33) 



Upon substituting these expressions into pT|) - p^ . pi 
(frai and (inil-dMl), followed by a transition to the hmit 
a/{u!eo) +00, it is verified that for a PEC surface 
only the terms lai in (P5)) ~ ([77|) are nonzero, for either 
polarization. 

Figure[2]compares the spatial correlation function (scf) 
of the normal component of the incident plus reflected 
electric field 



PEjkQAz;kQZQ) = 

(Ez{kQZo)E*{koZo 



fcoAz)) 



^{\Ez{koz„)\^){\Ez{koZo + koAzW 



(34) 



at koZo = tt/A for selected values of a/{Lueo) ^ 1 with 
the corresponding function for a PEC surface. Finite 
values of a/iujeo) are seen to cause the first zero crossing 
of '3i[pE^ (koAz)] to occur at smaller values of koAz com- 
pared to a PEC boundary. Scfs for Et and E (Figs. [3]and 
S]) show that corresponding differences for ^[pE^t^ {koAz)] 
between finitely conducting and PEC surfaces are less 
pronounced than for E^- Compared to a PEC surface, 
the damping of 5R[p£;j (/cqAz)] is qualitatively different 
from that for "SilpE^ikoAz)]. Also, finite conductivities 
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yield nonvanishing imaginary parts of the scf, indicating 
that E{koZo) and E*{koZo H- koAz) or their components 
are, on average, no longer in phase. This effect can be 
exploited as a means to measure surface conductivity. 
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FIG. 2: Scf of the incident plus reflected normal 
field Ez for selected values of a/{ujeo) at kozo = n/4 
as a function of separation fcoAz in normal direc- 
tion. Curves originating at ordinate value 1 repre- 
sent R[pEi,{koAz)]; curves originating at ordinate value 
represent Q[pE^{koAz)]. 
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FIG. 4: Scf of the incident plus reflected total (vec- 
tor) field E for selected values of a/{u!eo) at fcozo = 7r/4 
as a function of separation fcoAz in normal direc- 
tion. Curves originating at ordinate value 1 repre- 
sent 9?[pB(feoAz)]; curves originating at ordinate value 
represent Q[pE{koAz)]. 
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FIG. 3: Scf of the incident plus reflected tangential 
field Et for selected values of a/{u}eo) at kozo = 7r/4 
as a function of separation fcoAz in normal direc- 
tion. Curves originating at ordinate value 1 repre- 
sent 3?[pBt (fcoAz)]; curves originating at ordinate value 
represent Q[pEt{koAz)]. 
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B. Lossless isotropic dielectric medium 

As a second special case, we analyze the effect of the 
permittivity of a lossless isotropic dielectric medium on 
the scf. For brevity, we now limit the presentation to 
results for the amplitude of the vector field E only. 

For the refracted field. Fig. [S]shows that the permittiv- 
ity manifests itself by a decrease of the first zero crossing 
distance (correlation length) for 5R[p£;(fcoAz)], with asso- 
ciated shifts of the local maximum and minimum values 
toward lower values of kgAz. Also, the amplitudes of 
^[pEikoAz)] and 3[p£;(fcoAz)] increase with increasing 
e/eo- Similar findings apply to pe^ (koAz) and pEt (koAz) 
(not shown). The high sensitivity of psikoAz) to the 
value of e/eo, in combination with its insensitivity to 
koZQ, suggest that measurements of pEikgAz) may be 
used as a precision method for determining the refrac- 
tive index of a transparent substance. 
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FIG. 5: Scf of the refracted vector field E for selected 
values of ~ e/eo at arbitrary kozo as a function of 
separation koAz in normal direction. Curves originat- 
ing at ordinate value 1 represent ^[pE{koAz)]; curves 
originating at ordinate value represent '^[pE(koAz)]. 



For the incident plus reflected fleld, Fig. [5] shows 
psikoAz; koZo) in the half-space of incidence at kgZQ = 
7r/4, for selected values of e^ = e/eo. Comparing the 
asymptotic curve for e^ — > -f oo with Fig. |4] for a PEC 
surface, it is noticed that pE{koAz) is qualitatively sim- 
ilar, but quantitative differences exist, particularly for 
koAz < 1. Thus, the scf for the reflected field can be 
used to distinguish between conducting and high-k di- 
electric media, both of which exhibit high reflectivities 
making them otherwise difficult to discern in scalar local 
measurements. 
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FIG. 6: Scf of incident plus reflected vector field E for 
selected values of e^ ~ e/eo at koZQ = 7r/4 as a function 
of separation koAz. Curves originating at ordinate 
value 1 represent ^[pE{koAz)]; curves originating at 
ordinate value represent Q[pE{koAz)]. 



C. Directional incidence 

So far, the direction of incidence {9o, 0o) of the random 
field onto the interface was assumed to be uniform within 
the upper hemisphere {Qq = 27r sr, viz., —tt/2 < 9q < 
< </)o < tt). In practice, particularly in millimeter- 
wave and optical regimes, the wavevectors of the incident 
EM beams are often confined to be within a narrower 
solid angle 9o - A0o < 9o < Oo + A0o, 4>o - A<j)o < 4>o < 
(j)Q + A(j)o and scaling by 2A(j)Q/'K in azimuthal direction. 
For A0OjA(/)o 0, this approaches an unpolarized EM 
beam incident along the (0o, </'o)-direction. 

The spatial coherence along the refiected or refracted 
beam can be calculated as before, upon replacing the 
integrations du for the angular spectral averages by 

/cos(Aeo) -'^^'-'^S ^^'^ direction of specular reficction, 
— ^Oj the tangential and normal point separations are re- 
lated by Ax/Az = tan 6*0 = \J (1/u)^ — 1; along the di- 
rection of refraction, 9, their ratio is Ax/Az = — tanS = 
— Vl — I \/ {k/ko)'^ — 1 -|- M^. In general, narrow inci- 
dent and refiected beams fields do not interfere unless 
incidence is sufficiently close (with respect to the beam 
width) to the surface normal, whence the solid angles 
of the incident and reflected waves overlap partially or 
completely. 

To illustrate the effect of directionality of incidence 
on the scf, we consider an incident random field rep- 
resented by an angular spectrum of elevational width 
(fleld of view) 2A9q centered around the normal direc- 
tion (^0 = 0) with preservation ofthe azimuthal symme- 
try around this direction (A(^o = "') and random polar- 
ization (0 < ^/iQ < 2tt). This corresponds to incidence 
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from within a solid angle 27rA0o sr. For a PEC surface, 
Fig. [7| shows p£;(fcoAz) at selected values of A6o for 

-oo, 



^0^0 +00- It can be verified that for kozo 

PEikoAz;A9o^0) pE,{koAz; AOo - 

cos(fcoAz) 

On the other hand, for koZo — s- 0, 



0) 



PEikoAz;AdQ^O) pE^koAz; A9o ^ 0) 

sgn[cos(A:oAz)] 



(35) 



(36) 





FIG. 8: Scf of refracted vector field E at selected val- 
ues of 6*0 with e/eo = 2 as a function of separation koAz 
in normal direction. Curves originating at ordinate 
value 1 represent ^[pE{koAz)]; curves originating at 
ordinate value represent Q[pE{koAz)]. 



FIG. 7: Scf of incident plus reflected vector field E in 
normal direction in front of a PEC surface, at selected 
values of A^o for fco-^o — > +oo as a function of separation 

koAz. 

Fig. [5] shows corresponding results for refraction by an 
isotropic dielectric medium with e/eo = 2, demonstrat- 
ing qualitatively similar features for both 5R[p£;(fcoAz)] 
and ^[pE{koAz)], with pE^koAz) —+ exp(— jfco-^/e^Az) 
for A^o 0. In general, a medium with larger yields 
more rapidly modulated oscillations of its scf. 

Corresponding functions for the incident plus reflected 
field at kQZQ = tt/A are shown in Fig. [SI A general feature 
is that the oscillations of the scf become less regular when 
koZo decreases for a given value A^O: or vice versa; see, for 
example, the plot of pE{koAz) for A^o = 1 deg in Fig. [9l 
In the limit kozo 0, the scf tends again to the complex- 
harmonic square-wave function sgn[exp(— j'/cqAz)]. 
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FIG. 9: Scf of incident plus reflected vector field E at 
selected values of A9o with kozo = 7r/4 and e/eo ~ 2 
as a function of separation koAz in normal direc- 
tion. Curves originating at ordinate value 1 repre- 
sent ^[pE{koAz)]; curves originating at ordinate value 
represent '^[pE{koAz)]. 
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IV. CONCLUSION 

The reflection and transmission properties of an 
isotropic magneto-dielectric medium have been shown to 
have a significant effect on the spatial correlation func- 
tions of normal, tangential and total EM random vector 
fields, as a consequence of the EM boundary conditions. 
An analysis of the corresponding effects on the proba- 
bility distribution of the energy density of the total field 
for this configuration will be presented in a forthcoming 
paper. 

The effect of changes in the constitutive parameters on 
the correlation length is in general ambiguous, because a 
decreasing first-zero crossing distance of the scf is usually 
accompanied by an increase in its amplitude, \pE{kQAz)\, 
so that e.g. the two definitions investigated in Section V 



of [a| yield diverging tendencies for such changes. 

Since the above formulation is in terms of TE and TM 
refiection and transmission coefficients for plane waves, 
the analysis can be extended without effort to investigate 
the scf for more general stratified configurations with 
uncoupled eigenpolarizations of this kind, e.g., multi- 
layered stratified media as well as uniaxial anisotropic 
media, by substituting the appropriate functional forms 
of these coefficients. 
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